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We demonstrate that quantum well states (QWS) of thin Pb films are highly perturbed within the 
proximity of intrinsic film defects. Scanning Tunneling Spectroscopy (STM/STS) measurements 
indicate that the energy of these states have a strong distance dependence within 4 nm of the defect 
with the strongest energetic fluctuations equaling up to 100 meV. These localized perturbations show 
large spatially-dependent asymmetries in the LDOS around the defect site for each corresponding 
quantum well state. These energetic fluctuations can be described by a simple model which accounts 
for fluctuations in the confinement potential induced by topographic changes. 
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The curious yet prolific properties of thin Pb quan- 
tum films are rooted in its rich thickness-dependent 
quantum oscillation phenomena resulting from strong 
quantum size effects (QSE)^^. Robust quantum os- 
cillations of the underlying electronic states can have 
dramatic consequences related to growth, kinetics, and 
chemical reactivit}EM^. Most notably, this system has 
become a fruitful playground for investigating the in- 
terplay between quantum electron confinement and su- 
perconductivity. Initial reports demonstrated that the 
superconducting transition temperature, Tc, exhibits a 
pronounced thickness-dependent oscillation which is well 
correlated with the beating of the underlying quantum 
well states (QWS)PK^. 

While the initial focus was placed on the correlation 
between the superconducting gap and QWS, recent stud- 
ies have raised several questions^^. Surprising observa- 
tions of a pseudogap have brought into question the role 
of electron-phonon scattering as it relates to QSE^^. Fur- 
thermore, astonishing observations of superconductivity 
in 2 ML Pb films, which exhibit a strong Tc dependence 
on the underlying interface, independent of the QWS 
near Ef^ has also raised a question about a simple corre- 
lation between QWS and superconductivit}'^^. In order 
to understand how the QWS play a role in the above 
mentioned phenomena as well as other physical phenom- 
ena, a deeper and broader understanding of how QSE 
behave in these films is becoming a quandary that needs 
to be further addressee?^. 

We characterize commonly observed film defects 
within thin Pb films and demonstrate how the underly- 
ing quantum well states (QWS) can dramatically change 
within the vicinity of these defects. Using scanning tun- 
neling microscopy /spectroscopy (STM/STS), we show 
that the energetic location of all observed quantum well 
states can fluctuate as much as 100 meV depending on 
the lateral location relative to the defect. Within the 
vicinity of a defect, these modified states are highly 
spatially asymmetric resulting from a topographically- 
induced perturbation to the confinement potential. We 
describe these energetic fluctuations within a simple 
model where the QWS values are interpolated from dis- 



crete thickness-dependent values deduced for each film 
thickness, away from defects, measured by tunneling 
spectroscopy. 

All experiments were carried out using a home built 
STM operating at T = 6.5 K. Ptir tips and Ir tips were 
used for all experiments. STS was taken using a lock- 
in amplification scheme with a typical modulation fre- 
quency of 2 kHz and modulation voltages between 4 — 10 
mV. Pb was evaporated from a thermal evaporator and 
deposited onto either a clean Si(lll) — (7 x 7) surface or 
{V3 X V3)-Pb ("V3") interface held near liquid nitrogen 
temperatures^^. Globally flat Pb films are form ed after 
annealing cold films near room temperature^^^^. 

Pb films exhibit a rather high defect density when uti- 
lizing these growth conditions. Fig [TJa) illustrates an 
STM image of three different types of intrinsic defects on 
a Pb film (indicated by letters). All three defect types 
appear independent of film thickness (probed up to 25 
ML) as well as Si/Pb interface (7 x 7 or a/3). Here, we 
focus only on the properties of type A defects. 

Type A defects are oriented along one of three (Oil) di- 
rections. Each type A defect has a characteristic semicir- 
cular region of lower contrast, referred to as the "ditch", 
which is preferentially located along one side of the de- 




FIG. 1: (color online) (a) STM images of three types of defects 
(labeled by letter) seen on the surface of globally flat Pb films. 
White lines are used to help guide the eye. Vsampie = +1.5 V; 
It = 77 pA. (b)-(c) Atomic resolution of a type A defect indi- 
cating the ditch region (red arrow) and the spine region (blue 
arrow), (b) Vsampie = —300 mV; It = 25 pA (c) Laplacian 
flatting of an atomic image; Vsampie = —150 mV; It = 25 pA; 
(9 ML y/3 film). ). The number of atoms along the dashed 
line is indicated. 
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FIG. 2: (color online) (inset) Topography image of type A 
defects referencing spectroscopy taken in (a)-(b) Vsampie = 
+321 mV; h = 100 pA. (a) STS taken in a ditch region 
(blue) and in a spine region (orange) of an 18 ML Pb film 
(b) QWS energies derived from spatially resolved STS taken 
along a perpendicularly bisection line across the interfacial 
region between ditch and spine region of the defect indicated 
in (inset) by the arrow. 



feet for a given (Oil) orientation. These ditches have a 
typical corrugation of 0.5 — 0.9 A and radius of 2.5 — 4 
nm. For a given (Oil) defect orientation, there is an 
abrupt interface between the ditch and a localized region 
of higher contrast, which we refer to as the "spine," that 
extends along the dislocation line. Both spine and ditch 
regions maintain similar contrast at all observed imaging 
conditions regardless of which part of the unit cell is im- 
aged. Therefore, the resultant contrast is topographic in 
origin. 

Figjljb) illustrates atomically resolved images of a type 
A defect. A clear 1x1 pattern corresponding to an 
fee close-packed surface with a periodicity of 0.35 nm 
can be seen. As shown in Figjljc), a Burger's circuit 
taken around the defect exhibits one extra atom along the 
spine side of the defect corresponding to a Burgers vector 
b = I [Oil]. While this is equivalent to a single perfect 
dislocation, a dislocation of this type cannot account for 
experimental observations. According to Frank's Energy 
Criterion, a perfect dislocation with bi will dissociate 
into partial dislocations b2 and hs if {b1 > ^2 + ^3)^^- 
a perfect dislocation of the given form, this results in two 
Shockley partial dislocations of the form b^ = |(112). 
This dissociation produces a stacking fault on the surface 
plane, but will not produce a spine or a ditch. A possi- 
bility which will produce a ditch-like feature is a stacking 
fault which is situated along a close-packed plane off-axis 
(e.g. (Ill)) from the surface normal (e.g. 70.5°) which 
extends into the surface. The displaced atoms beneath 
the surface can produce an overall contrast change which 
would result in the appearance of a ditch but maintain 
atomic order as viewed from the surface. We believe 
this is responsible for the ditch feature. However, favor- 
able dissociations of a perfect dislocation along an off- 
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FIG. 3: (color online) (a) - (c) Relative shifts in QWS energies 
derived from STS for 3 separate defects (symbols) . refers 
to the energy location of a QWS subtracted from the state 
energy away from the defect. Each color refers to a specific 
defect. Solid lines refer to the calculated energy shifts from 
the interpolation model, (d) Reference line profiles for the 
three defects used to calculate the energy shifts in (a)-(c) 
(reference indicated in Fig 2(a)). 



axis close packed plane or a Frank partial (bp = |[111]) 
would alone not conserve the measured Burgers vector. 

Combinations of dislocations can dissociate and inter- 
act t hereb y producing features similar to both ditch and 
spin d^-^1^^ 4 For the PtNi(lll) surface, similar ditches 
have been attributed to dissociations of a subsurface 
dislocation network which produces an off-axis stacking 
fault along a close-packed plane which is bounded by 
other partial dislocations. STM images of these types 
of ditches are very similar to what is observed for a type 
A defect. Furthermore, it was demonstrated that the 
depth of the dislocation network is related to the overall 
depth of the ditch and shape of the ditch as seen with 
STM. This can explain the variation in the corrugation of 
our ditch. However, the exact structure which produces 
these PtNi ditches are not accompanied by a vacancy row 
(spine) and the ditches are not hemispherically shaped. 
Nevertheless, this again suggests that multiple disloca- 
tions, with subsurface components, are involved in this 
structure. 

Other possible dislocations where two perfect dis- 
locations along different close-packed planes dissoci- 
ate and interact to form a sessile dislocation (stair- 
rod dislocation/Lomer-Cottrell lock) are intriguing 
option^. These structures have subtle differences with 
what is observed in our case. Without further experi- 
mental data or simulations which can produce informa- 
tion about participating subsurface dislocations, it is not 
directly evident what exact dislocation(s)/dissociations 
beneath the surface produce this defect. 

Scanning tunneling spectroscopy (STS) within the 
vicinity of these defects reveals well-defined states re- 
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FIG. 4: (color online) Topography and constant-current dl/dV images {Vmod — 4 mV). (a) Topography Vsampie — +470 mV; 
It = 95 pA. (b) Constant-current dl/dV image Vsampie = +321 mV; It = 100 pA, (scale: 0.8 — 2.2 V) (c) Constant-current 
dl/dV image Vsampie = +470 mV; It = 95 pA, (scale: — 1.4 V). Arrows indicate specific defects which show different 
characteristic spatial distributions. All Pb films were 18 ML thick grown on a Pb-Si interface. 



suiting from quantum confinement^ However, there 
is a vast variation in QWS energies within the vicin- 
ity of a defect depending on the probed location (Fig 
[2ja)). QWS exhibit lower overall energies when probed 
in the spine (orange) and higher overall energies when 
probed in the ditch (blue). There is a clear spatially- 
dependent trend in the energetic shifts of the QWS (Fig 
[2|b)) which shows the strongest modulations near the 
interfacial region. Each color/symbol corresponds to a 
particular state for the given film thickness. Less than 4 
nm from the interface, specifically within the ditch region 
of the defect, there is a monotonic increase in all state 
energies as the orthogonal distance decreases between the 
probed region and the interface. At the interface, there 
is an abrupt decrease in all the state energies which is far 
below the average energy value for each state for a given 
thickness away from the defect. These values monotoni- 
cally increase as the orthogonal distance from the spine 
region increases until all state energies are restored to 
their typical value away from the defect. All energetic 
shifts are localized within a radius of ~ 4 nm from the 
center of the interfacial region. 

Fig|3|a)-(c) illustrate the measured energy fluctuations 
across three different type A defects on a 18 ML film. 
Here AE represents the state energy at the given loca- 
tion subtracted from the average state value of the film 
away from the defect. Negative distance values reference 
the ditch region (Figjsj^d)). Energetic shifts within the 
ditch region are persistently larger in magnitude, for each 
state, as compared to the spine region. Furthermore, the 
shifting behavior shows a weak state dependence, espe- 
cially for the ditch region, where the overall shifting at 
a given point increases with absolute energy. All type 
A defects show the same energetic shifting behavior in- 
dependent of the defect orientation and underlying in- 
terface. All spectroscopic points for a given defect were 
taken using similar stabilization conditions in order to 
avoid Stark shift effects created by the tipP^. 

The overall change in state energy can be understood 



within the context of a simple perturbation to the con- 
finement potential. By interpolating the l-th QWS Ei^n 
energy value for a given discrete thickness n, dEi^n/dz 
can be extracted for non-integer thickness values where 
z represents any real- valued thickness^. The simulated 
changes in state energies can be calculated by multiply- 
ing the deviations in topography near the defect, Az, by 
dEi^n/d^^' The solid lines in Fig[3|a)-(c) represent the 
calculated changes in energy for each defect assuming the 
absolute height away from the defect is the zero reference 
(Fig|3|d)). This simple model reproduces the observed 
experimental trends such as the state dependence of the 
energetic shifting as well as the overall order of magni- 
tude of the changes in state energy in the defect regions. 

The overall enhancement of the LDOS near a defect 
site shows a highly non-homogenous spatial distribu- 
tion at particular energies. Fig [4] illustrates constant- 
current dl/dV images and the correlated topography im- 
age taken below {E +321 meV; AE ^ -60 meV) and 
above {E = +470 meV; AE ^ +90 meV) (Fig a given 
QWS for a typical 18 ML ^/3 film. Strong enhancement 
of the LDOS can be seen near type A defects both be- 
low and above the given state energy {Eqws = +379 
meV). For downward energy shifts, the enhanced LDOS 
in the spine region shows a more isotropic and localized 
distribution. The enhancement spans the length of the 
interfacial region and is localized ~ 2 nm within the spine 
region. For the ditch region, the LDOS exhibits similar 
enhancements. However, unlike the spine region, there 
are pronounced geometrical features with strong intensi- 
ties in the ditch region of each defect (indicated by ar- 
rows). The number of lobes and the overall shape of these 
features varies depending on the diameter and apparent 
depth of the ditch. However, the geometrical shape for 
a given defect is nearly identical above all probed states. 
This is unlike the case where localized lobe features, at- 
tributed to quantum confinement, change as a function 
of energ^p^^. These lobes have a characteristic diameter 
between 1.5 — 4 nm. There is some indication that these 
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lobes may be quantized in multiples of Ai?/2. We believe 
these spatial distributions may be understood within the 
context of electron scattering within proximity of the dis- 
locations of the defect. However, this requires a more 
detailed model than the one presented. 

A description of the resultant dipole moment em- 
anating from a defect site was first predicted by 
Smoluchowski^^. Within the context of the jellium 
model, a dipole moment can result from an abrupt change 
in morphology which leads to a smearing in the charge 
density on the order of Ai?/2. However, a more com- 
plete picture of how this resultant polarization is pro- 
duced requires a detailed account of the modified bond- 
ing of both atoms and vacancies within the vicinity of the 
defect^^. Few measurements have been able to detect this 
polarization^^ For these cases, the measured fields 
can be described by a step-edge model. Descriptions of 
more complex defects have eluded theoretical modeling 
because of the intensive resources required. In context of 
this work, an overall modulation of the QWS naturally 
results in a modulation of the work function which should 
produce a local polarization. The field should be local- 
ized on the order of Ai? /2 within the defect along the off- 
axis close packed plane of the defect where Ai?/2 ^ 0.5 



nnPH. This local field may further perturb the QWS 
within the vicinity of the defect. This may explain why 
the downward shifts in the spine region, within 1 nm 
of the interfacial region, are much larger than predicted 
from our model. Nevertheless, from spectroscopy alone, 
we cannot distinguish between QWS fluctuations result- 
ing from topographical perturbations and those produced 
by dipolar contributions. 

In conclusion, we demonstrate that within the vicinity 
of intrinsic defects there exists large modulations in the 
QWS {AE ^ 100 meV). These fluctuations can be de- 
scribed by a topographically-induced perturbation to the 
local confinement potential. LDOS maps reveal that at 
the corresponding local QWS energy near the defect sites 
that the LDOS has a rather complex distribution. These 
large fluctuations in the overall LDOS should be consid- 
ered as these Pb systems become increasingly utilized for 
a variety of experiments involving kinetics, surface chem- 
istry, and superconductivity. 
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